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Abstract 

We consider a mixed boundary-value problem for the Poisson equation in a plane 
(^ [ thick junction fig which is the union of a domain 0,q and a large number of e— periodically 

^^ ' situated thin rods. The nonuniform Signorini conditions are given on the vertical sides of 

the thin rods. The asymptotic analysis of this problem is made as e — > 0, i.e., when the 
number of the thin rods infinitely increases and their thickness tends to zero. With the 
help of the integral identity method we prove the convergence theorem and show that the 
nonuniform Signorini conditions are transformed (as e — > 0) in the limiting variational 
inequalities in the region that is filled up by the thin rods in the limit passage. The 
-^ ^ existence and uniqueness of the solution to this non-standard limit problem is established. 

^ • The convergence of the energy integrals is proved as well. 

^' 

Ci . Key words: homogenization; thick junction; Signorini boundary conditions; variational 

inequalities 
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1 Introduction and statement of the problem 

(N 

C^ ■ Boundary-value problems in thick junctions are mathematical models of widely used engineer- 

ing and industrial constructions as well as many other physical and biological systems with 

O ■ very distinct characteristic scales. In recent years, a rich collection of new results on asymptotic 

analysis of boundary- value problems in thick multi-structures is appeared (see [lj-[8j). 

In this paper we homogenize the Signorini problem in a thick plane junction of type 2:1:1 



>< 



^ I using the integral identity method developed in [9l |T0] . 

A thick junction (or thick multi-structure) of type k : p : d is the union of some domain 
in M", which is called the junction's body, and a large number of e-periodically situated thin 
domains along some manifolds on the boundary of the junction's body (see Fig. 11.11) . This 
manifold is called the joint zone. Here £ is a small parameter, which characterizes the distance 
between neighboring thin domains and their thickness. The type k : p : d oi a thick junction 
refers respectively to the limiting dimensions of the body, the joint zone, and each of the 
attached thin domains. 

This classification of thick junctions was given in [II]-[T6] and [HI HO], where rigorous 
mathematical methods were developed (homogenization, approximation, asymptotic expan- 
sions) for analyzing the main boundary-value problems in thick junctions of different types. 
It was pointed out that qualitative properties of solutions essentially depend on the junction 



type and on the conditions given on the boundaries of the attached thin domains. In addition, 
as it was shown in |17j such problems lose the coercitivity as e — >■ and this creates special 
difficulties in the asymptotic investigation. It should be noted that papers [IHl IS] are the first 
papers in this direction. 

For the first time a problem known now as the Signorini problem was considered by Sig- 
norini himself in [20]. The sense of the Signorini boundary condition consists in a priori 
ignorance which of the boundary conditions (Dirichlet or Neumann) are satisfied and where. 
Many interesting problems in applied mathematics involve the Signorini boundary conditions. 
Applications arise in groundwater hydrology, in plasticity theory, in crack theory, in optimal 
control problems, etc. (see [21j). Such of these problems as can be recast as variational inequal- 
ities become relatively easy to study (see [21]- [23]). Asymptotic investigations of variational 
inequalities in perforated domains were made in [21 

The results of this preprint will be published in 



1.1 Statement of the problem 

Let a, I be positive numbers, h be a fixed number from the interval (0, 1), and A^ be a large 
positive integer. Define a small parameter e = -^. A model plane thick junction ^2^ (see Fig. II. 11) 
consists of the junction's body 

VLq = \^x = (xi,X2) G M^ : < Xi < a, < X2 < 7(3^1)} , 

where 7 G C^([0, a]), and a large number of the thin rods 



G,{e) 



X : 



T-^' + r 



<|, x2GH,0][>, i = 0,l,...,A^-l, 



AT-l 



i. e. n, = noU Ge, where G^ = IJ Gj{e). 

j=0 




B,-(0,a)x(-l,0) 



Figure 1: A model plane thick junction fig of type 2:1:1 



The discrete parameter e characterizes the distance between the rods and their thickness 
that is equal to eh. Obviously, the thin rods fill out the rectangle Dq = (0, a) x (— /,0) in the 
limit passage as A^ — >• +00 {e ^ 0). 

Let us denote the union of vertical sides of the thin rods Ge by Se] the union of bases of 
the thin rods will be denoted by Fg. 

In Q^ we consider the following boundary value problem 

-AUe{x) = f{x), X e fie, 

Ue{.x) < g{x), duUe{x) < ed{x), X G Se, 

(1) 

{us{x) - g{x)) {dyUe{x) - ed{x)) =0, x E S^, 

^ ue{x) = 0, xe r„ d^usix) = 0, xedQs\ {Se u r,), 

where d^ = d/dv is the outward normal derivative, /, g, d are given functions. 

We assume that / G L^(fii), where fii = fio U Dq, the function d belongs to the Sobolev 
space H^{Dq), and 

g G H\Do; h U h) = [v e H\D^) : v\,^^,^ = O} , 

where Ii = {x : Xi G (0, a), X2 = —1} , Iq = {x : Xi G (0, a), X2 = 0} . 

Our goal is to study the asymptotic behavior of the solution Ue to problem ([T]) as e — > 0, i.e., 
when the number of the thin rods infinitely increases and their thickness tends to zero. 

2 Definitions of the weak solution and its existence 

In the Sobolev space H^{VLe] T^) = {m G H^{Vte) '■ ^lr = O} , we define subset 

Ke = {ip e H^{Vte]Te) : ip\g^ < g\g^ a.e. on S'J, 

where ip\g denotes the trace of a Sobolev function ip on a. curve S. Obviously, K^ is a closed 
and convex set for every fixed value of e. 

Let us suppose the existence of a classical solution to problem ([1]). We can regard that 
g = in Qq. Multiplying the equation of problem ([T]) by the function {u^ — g), integrating by 
parts in fl^ and taking into account the boundary conditions for u^, we obtain 

/ \/ue-Vu£dx+ / Vus-V{u£ — g) dx = / fuedx+ / f{ue — g)dx + e / d{x){u£ — g) ds, (2) 

where Vw ■ Vw = f] ^§^- 
j=i ^ ■' 

Now we take any function tp & K^ and multiply the equation of the problem ([T]) by {ip — g). 

Similar as before we get 

/ Vwe ■ Vv? dx + Vug ■ V{(p — g) dx = 



fipdx+ f{ip — g)dx + ed{x){ip — g)ds+ {d,yUs — ed{x)){ip — g) ds. (3) 

Ho Ge Se Se 

Since d^Ue{x) < ed{x) and ip{x) < g{x) a.e. in S^, 

{d^u,-ed{x)){ip- g)ds>0. (4) 

Taking into account (JH), it follows from equality (E]) that 

/ Vwe ■ V(p dx + Vue ■ V(v9 — g) dx > f(pdx+ / f{(p — g) dx + e / d{x) {(p — g) ds (5) 

Q,0 Ge CIq Ge Se 

Definition 1. A function Ue € Kf, is called a weak solution to problem (Qp if it satisfies the 
integral equality ^ and integral inequality ^ for arbitrary function if G K^. 

Another definition is as follows. 

Definition 2. A function m^ G K^ is called a weak solution to problem (J\) if it satisfies the 
integral inequality 



VUs ■ V((y9 — Ue) dx > f{(p — Us) dx + E / d{x){(p — Ue) ds V (yS G K^. (6) 

Let us show that these definitions are equivalent. Subtracting equality ^ from inequality 



), we arrive at (jH]). Setting ^p = I ' „ ' in ([H]), we have 

- / \Vue\^ dx+ / 'S/Ue-'S/{g — Ue)dx > — I fuedx+ I f{g — Ue)dx + e I d{x){g — Ue) ds. (7) 



Ho Ge no Ge 



Putting (^=<^' „'in([6]),we get the reversed inequality 

iWuef dx + jVue ■'V{us — g) dx > fuedx+ f{us — g) dx + e d{x){ue — g) ds. (8) 



Qq Ge fio Ge 



This means that (I2]) holds. Setting o? = s , , ^' ^ ' in ([6]), where ib is arbitrary 

t ^ + ^£-5', a; G Ge, 

function from Ks, we get ([5]). 

It is well known (see for instance [2T]-[23]) that there exists a unique solution of inequal- 
ity dn]) for any fixed value of e. 



3 Auxiliary uniform estimates 

To homogenize boundary-value problems in thick multi-structures with the nonhomogeneous 
Neumann or Fourier conditions on the boundaries of the thin attached domains the method of 
special integral identities was proposed in [9l [10] • For our problem this identity is as follows 
(see [TtJi Lemma 1]) 



-— / f dx2 = V dx — e Y i — j dx-^v dx W v E H^ [Ge) , 



(9) 



where Y{^) = —C, + [^] + |, [^] is the integral part of ^. 

Using ([9]) and taking into account that maxjR |y| < 1, we get 

\\v\\ms.) < Cis--^ \\v\\hhg.) yve H\G,). (10) 

Remark 1. Here and in what follows all constants {Ci} and {cj} in inequalities are indepen- 
dent of the parameter e. 

Also it will be very important for our research the following uniform estimates. 

Lemma 1 ([14j). The usual norm ||'u||/^i(q^) = ( J^ {\V u\'^ + v?) dx\ in H^iVt^^Ve) and a 
norm || ■ ||e, which is generated by scalar product 



{U,v)e 



/ Vu-Vvdx, u,v e H^i^s^Te 



are uniformly equivalent, i.e., there exist constants Ci > and ^q > such for all e G (0,£:o) 
and for all u G if ^(fi^, F^) the following estimates 

\\u\\e < \\u\\m(^a,) < Ci\\u\\^ (11) 

hold. 

Remark 2. In fact, in LemmaUl the following Friedrich inequality 

IklU^cn.) < C2II VM|U2(n,) Vug H\Q„ F,) (12) 

was proved. 

Using the Cauchy-Bunyakovsky integral inequality and Cauchy's inequality with 6 > 
{2ab < 8a? -\- 5~^6^ for any positive numbers a and 6), with the help of (ITUl) and (IT^ we deduce 
from (E]) that 

/ IVWeP dx < Co{6i +62 + 63) II VMe||^2(n,) + 

+ Ci(l + 6^') Ml^^no) + C2(l + '^2"') ll/lli^(n.) + C3(l + S3') Ml.^^^y (13) 

Choosing 61, 62, 83 so that Ci{6i + 62 + 63) < |, we have 

j \Vu,\'dx < C4 (iml^^m) + MmiDo) + MmiDo)) ■ (14) 

By virtue of flTTj) we obtain from flT^ the following uniform estimate 

WueWn^n,) < C3. (15) 



4 Convergence theorem 



In the sequel, u denotes the zero extension of a function u defined on Gg to the rectangle Dq, 
which is filled up by the thin rods in the limit passage as e — i> 0, namely: 



uix] 



u{x), X E Gs, 

0, XEDo\ Ge. 



If M G H^{Gs, Fg), then due to the rectilinearity of the boundaries of the thin rods u belongs 
to the anisotropic Sobolev space 

H^^\Do]Ii) = {v e L^{Do) : 3 weak derivative d^^v e L'^{Dq) and v\j^ = O} (16) 

and d^^u = 8x2^ a-e. in Dq. 

Theorem 1. The solution Ue to problem ([2]j satisfies the relations 

Up\r,.. -^ Uq weakly in H^{VLq), 

'-^ Huq weakly in W^'^{Dq,Ii)^ )• as e -^ 0, 






(17) 



and the function Uo{x) 



->■ weakly in LF'{Dq) 

1' ' is a unique solution of the following problem 



-A<(x) 
-hdl^x2^Q{x) 
Uo{x) 



a;2a;2 



d^u^ix) 

Uq{xi,-1) 

Uq{xi,0) 



which is called the homogenized problem for ^\). 

Furthermore, the following energy convergence holds 



= fix), 


X G rio. 


<hf{x) + 2d{x), 


X G Z^o, 


< 9ix), 


X G Z^o, 


2f{x) + 2d{x)) = 0, 


X G Do, 


= 0, 


X edQo\ Jo, 


= 0, 


Xi G [0,a], 


= u^{xi,0), 


Xi G [0,a], 


= hdx^Uoixi^O), 


Xi G [0,a], 



(18) 



\im E,{us) = Eo{uo) 



where 



Esiug) = / \\/Us\'^ dx, Eo{uo) = / \\/uq\'^ dx + h / \dx2UQ\'^ dx. 
Before the proof of Theorem [1] we investigate the homogenized problem (fTSll . 



4.1 The solvability of the homogenized problem 

We see that the homogenized problem (1181) is a non-standard boundary-value problem that 
consists of the Poisson equation in the junction body Qq, the variational inequalities in Dq and 
the transmission conditions in the joint zone Iq. Therefore, at first we give the definition of 
the weak solution to this problem and then with the help of general approach in the theory of 
variational inequalities we prove the existence and uniqueness. 
Let us introduce partly anisotropic Sobolev space 

n{Qu Do; Ii) = {ue L\Q,) : 3 d,,u G L\Q,), u\^^ e H\Qo), m|^^, G H''\Do; h)} , 

where H^'^{Dq]Ii) is defined in (iT6ll . It follows from properties of anisotropic Sobolev spaces 
(see [3n]) that the traces of the restrictions u"*" := u\^ and u~ := m|^ on Jq are equal. In 
addition, since traces of functions from 7i{Qi, Dq] Ii) vanish on Ii, there exists a constant Cq 
such that 

f u'^dx<Co(f \Vu+\'^dx+ f \d.^^u~\^dx) yuenini,Do]ii). 
Jui ^Jno Jdo ' 

In 7^(fii, -Do; ^i) we introduce a norm || ■ Ht^, which is generated by a scalar product 
{u,v)n= / Vu'^ ■ Wv^ dx + h / dx2W dx^v' dx, u, v E H{^i, Dq; Ii). 

We now define subset Kq = {ip E TiiVti, Dq; Ii) : (p~ < g a.e. in Dq} in Ti.{Qi,DQ; Ii). 
Obviously, Kq is a closed and convex set. 

Definition 3. A function uq G Kq is called a weak solution of problem fT^) if it satisfies the 
integral equality 

/ Vuq -Vuq dx + h / dx2UQ dx2{uQ — g) dx = / fu^ dx + h / f {uq —g) dx + 2 / d{uQ —g)dx, 

^0 -Do i^o Do -Do 

(19) 
and the integral inequality 

I Vuq -Vipdx + h I dx2UQ dx^i'-P — g) dx > / f(pdx + h / f {(p — g) dx + 2 / d{(p — g)dx (20) 

f^o Do no Do Do 

for arbitrary function if G Kq. 

If there exists a classical solution to the homogenized problem (llSp . then relations (I19p and 
(I^Uj) can be obtained by the same way as relations ([2]) and i^ in Definition [H Similarly as we 
proved the equivalence of Definitions [1] and [2], we can show the equivalence of Definition [3] to 
the following definition. 

Definition 4. A function uq G Kq is called a weak solution to problem fT^) if it satisfies the 
integral inequality 

/ Vuq ■ V{(p — Uq) dx + h dx2UQ dx2{'^ — Uq) dx > 
Jno J Do 



> / f{ip — UQ)dx + h / f {if — Uq) dx + 2 d{ip — Uq) dx (21) 

Jqo J Do J Do 

for arbitrary function ip G Kq. 



We now give the third definition of the weak solution to problem (fTS!) . 



Definition 5. A function uq G Kq is called a weak solution to problem fl^) if it satisfies the 
integral inequality 

/ Vifi ■ V{(p -u^)dx + h / d^^V dxziv - ^o") dx > 

JUo J Do 



> / f{ip — Uq) dx + h f {ip — Uq) dx + 2 d{ip — Uq) dx (22) 

Juo J Do J Do 

for arbitrary function ^p G K^. 

Let us prove that Definition H] and Definition [5] are equivalent. Adding the inequality 

/ V{ip~u^) ■V{ip-u^)dx + h dx2i^~UQ)dx2iv-UQ) >0 {(p e Kq) 

Juo Jdo 

to inequality (I7T1) . we get (1^ . Now we take any ip G Kq. Setting ip = uq + t{%lj — Uq) G Kq 
(for any t G [0, 1]) in inequahty (l22il . we obtain 

/ V iu'^ + t{%l) - u'^)) ■ V {%l) - u^) dx + h j dx.,{uQ + t{ij - Uq)) dx.,{ij - Uq) dx > 
Jno J Do 



> f{^-u^)dx + h f{^-UQ)dx + 2 d{i)- Uq) dx. (23) 

Jq.0 J Do J Do 

Passing to the limit in (J23ll as t ^ 0, we arrive at ( |2T1) . Thus, all Definitions [3], H] and [5] are 
equivalent. 

We can re-write inequality (I2T1) in the following form 



(m, ip-u)n>{F,ip-u) V ^ G Kq, (24) 

where F is a linear continuous functional on 7i(r2i, Dq] Ii) defined by the formulae 

{F,w) = / f w^ dx + h f w^ dx + 2 dw~ dx for all w G 7Y(i7i, Do! -^0- 

Jno J Do J Do 

Using the theory of variational inequalities in Hilbert spaces (see [211 Sec. 2]), we can state 
that there exists a unique solution of the inequahty ( 12^ and consequently of the homogenized 
problem (IT8|) . 

4.2 The proof of Theorem [1] 

1. From (dS]) it follows that \\ue\\m(no) < C's, WusWl-^^Do) < C's and Wd^MeWv^^Do) < C's, 
i = 1, 2. Therefore we can choose a subsequence {e'} C {e} (again denoted by e), such that 

"el no ~^ ^t weakly in H^{9-q), 

Us -^ huQ weakly in L?'{Dq), ^ as e -^ 0, (25) 

dxiUe -^ 7i weakly in L^{Do), i = 1,2, 

8 



where Uq, Uq, 71, 72 are some functions which will be determined later. 

At first we determine 72. Take any function ip G C^{Dq) and perform the following calcu- 
lations 



dx2U£ Tp dx = I dx2U£ i' dx = I S^jMe Tp dx = — j Us 8x2^' dx = — Us dx2''P dx. 

Do -Do Ge Ge Do 

Passing to the limit in this identity, as £ ^ 0, we obtain 

/ ^2^dx = -h I u^ dx2ilJ dx, V^ G C^{Dq), (26) 

Jdo Jdo 

whence it follows that there exists a weak derivative dx2U^ and 72 = hdx2UQ a. e. in Dq. 
Now let us find 71. Consider the function 

[eYi{f}ip + g, xeGs, 
where ^1(0 = ~^ + [^]- K is easy to see that $ G -ft'e and 

V($ -g)= (-^ + eY, (^) dx,^, eYi (^) dx2^) , x e G,. 
Substituting the function ^ — g into the integral inequality ([5]) for solution Ue, we get 

[ (-dx,Ueip + eYi (—)vue-Vip) dx>e f Yi (—) f ip dx - ^ ^^^ ' f dipdx2, 

G. G, 5± 

where the sings " + " or " — " in Sf" indicate the union of the right or left sides of the thin rods 
respectively. With the help of flTOj) and flTSj) we deduce from previous inequality the estimate 

f d^.ijdx <en \Y, (^) (Vu, -Vij-f ij)\dx + ^^^ ^^^ f \d^\dx2^ < 

<eCi {\\VUe\\L^G,)\\^HLHG,) + \\f\\LHG,)MLHG,) + 4d\\L^{S.)ML^(S,)) < 
<£Ci {\K\\Hi(n,)\\i'\\HHDo) + \\f\\LHn^)M\LHDo) + \\d\\HHDo)MH^iDo)) < £C2, 

from which, passing to the limit as e ^ 0, we get J^ ■jitp dx = for all ip G C^{Dq), V' > 0. 
This means that 71 = a. e. in Dq. 

2. Let us show that the traces of the functions Uq and Uq on /q are equal. By virtue of the 
compactness of the trace operator and the first relation in (1251) . we have 

Us{xi,0) -^ Uq{xi,0) inL^(0,a) as e —> 0. (27) 

Consider the following equality 

Ueixi, 0) = Xh i—j Ueixi, 0) for a. e. xi G (0, a), (28) 

9 



where X/i(0) ^ G IR, is the 1-periodic function, defined on the segment [0, 1] as follows: 

fl |£_ 11 < ^ 

"'•©-jo. |<|{-i|<l. 

It is known that Xh {—) — ^ h weakly in L^(0, 1) as e ^ 0. Using this fact and flTTj) . we 
obtain that the right-hand side in ( !28l) converges to h Uq{xi, 0) weakly in L^(0, a). On the other 
hand, 

a 

1 /" _ If ' 

U!r{xi,0)ip{xi)dxi = - / Us{x)4'{xi)dx+- j {x2 + l)dx2Ue-'^{xi)dx V'?/' G C^(0, a). (29) 

Do Do 

Passing to the limit in fl^Ul) as £ ^ and taking (1261) into account, we have 



h / Uq{-,Q)iI){xi) dxi = - Uqi/j^xi) dx + - {x2 + I) O^^Uq ip{xi) dx W i/j e C^(0,a), 

Do Do 

whence it appears 

I pa 

u+i;0)iP{x,)dxi= u^i;0)^{xi)dx, V V^ G Co^(0, a), 
Jo 

i.e., Uq{xi, 0) = Uq{xi, 0) for a.e. xi G (0, a). 

Similarly we can prove that the trace Uq\i^ is equal to zero. 
Thus, the results obtained above mean that the function 

( u^, X e Qq, 

belongs to the space Ti.{Qi, Dq; Ii). 
3. Let us add the inequality 

/ V{ip-Ue) ■V{ip-Us)dx+ dr^^iv - Us) d^^i^ - Us) dx + dr,^Ued^^Usdx > 0, 

J S7q J G^ J G^ 

where ^p is arbitrary function from C^(fii) such that Lp\i =0 and '^ < g \\y Dq (obviously 
f\n^^ Ks), to inequality (E]). We get 

/ V(p-V{(p-Us)dx+ d^^Us d^^ip dx + d^^ip d^^i^p - Ue) dx > 

•/iiQ JG^ JG^ 

> / f{(p — Us)dx + e / d{x){(p — Ue) ds, 



10 



which we can re-write with the help of (Q in the following view 

/ V(p ■ V{ip - Ue) dx+ d^^Us d^^ifdx + Xh ( — ) d^^Lpd^^if dx - d^^^p d^^u^ dx > 

Ho Do Do Do 

> f {<^-Ue)dx+ Xh y—J ffdx- fue dx+ 

no Do Do 

+ T Xh( — )d(pdx-- dusdx-— Y ( — ) d^^{d{ip - u^)) dx. (30) 



Do Do 



Passing to the limit in (1301) as £ — * and taking into account results obtained above, we 
obtain the following integral inequality 

/ Vy? ■ V(v9 — Uq) dx + h dx2^dx2{'~P " '^o) dx > 

Ho Do 

> / f{(p — Uq) dx + h / f {(p — Uq) dx + 2 I d{(p — Uq) dx (31) 

!^o Do Do 

for any function (p E Ki = [ip E C^{ili) : ip \j = 0, ip < g in Dq^ . 

Since the set Ki is dense in Kq, the integral inequality (13T|) holds for any function ip E Kq. 
This means that the function uq is the unique solution of inequality (|2T1) (see Definition H]) and 
also it is the weak solution to the homogenized problem flTHl) . 

Due to the uniqueness of the solution to problem flTSj) . the above argumentations are true 
for any subsequence of {e} chosen at the beginning of the proof. Thus the limits ( TT71) hold. 



4. /^From equalities ([2]) and (1191) it follows that 
Eeiue) = I \S/Ue\^ dx = I Vus-Vgdx+ I fui;dx+ / f{ue~g)dx + e I d{x){ui; — g)ds^ (32) 

Qg Gs ^0 G^ Se 

Eoiuo) = / {VuqI"^ dx + h / \dx2UQ l"^ dx = h / dx2Uo 9x2 9 dx+ 

iio Do Do 

+ / fuQ dx + h f {uq — g) dx + 2 / d {uq — g) dx. (33) 

no Do Do 

Passing to the limit in (I32p similarly as we made this in (I30p and taking into account (I33p . we 
obtain lim^^Q E^{us) = Eq{uq). The theorem is proved. 
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